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(Dated: May 12, 2005) 

We present a measurement of the fraction /+ of right-handed W bosons produced in top quark 
decays, based on a candidate sample of ti events in the lepton+jets decay mode. These data 
correspond to an integrated luminosity of 230 pb _1 , collected by the D0 detector at the Fermilab 
Tevatron pp Collider at y/s = 1.96 TeV. We use a constrained fit to reconstruct the kinematics of 
the ti and decay products, which allows for the measurement of the leptonic decay angle 0* for each 
event. By comparing the cos#* distribution from the data with those for the expected background 
and signal for various values of /+, we find /+ = 0.00 ± 0.13(stat) ± 0.07(syst). This measurement 
is consistent with the standard model prediction of /+ = 3.6 x 10~ 4 . 



PACS numbers: 14.65.Ha, 14.70.Fm, 12.15.Ji, 12.38.Qk, 13. 38. Be, 13.88.+e 



The top quark is by far the heaviest of the known fermions and is the only one that has a Yukawa coupling 
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of order unity to the Higgs boson in the standard model. 
The top quark is also unique in that it decays through 
the electroweak interaction before it can hadronize. In 
the standard model, the top quark decays via the V — A 
charged current interaction, and almost always to a W 
boson and b quark. We search for evidence of new physics 
in the t — > Wb decay by measuring the helicity of the W 
boson. The W bosons produced from these decays are 
predominantly in either a longitudinal or a left-handed 
helicity state with fractions fo and /_, respectively. For 
any linear combination of V and A currents at the tWb 
vertex 0, 



fo 



2M 2 , 



0.703 ±0.012 



(1) 



where m t is the mass of the top quark for which we use 
175 ±5 GeV (consistent with the world average @|), Mw 
is the mass of the W boson, and rrib is the mass of the bot- 
tom quark. In this analysis, we fix fo at 0.7 and measure 
the positive helicity (or right handed) fraction / + . In the 
standard model, /+ is suppressed by a factor of (mf,/m t ) 2 
and is predicted at next-to-leading order to be 3.6 x 10 -4 
A measurement of /+ that differs significantly from 
this value would be an unambiguous indication of new 
physics. For example, an / + value of 0.3 would indicate 
a purely V + A charged current interaction. A possible 
theoretical model that includes a V + A contribution at 
the tWb vertex is an SU{2) L x SU{2) R x U Y (l) exten- 
sion of the standard model Direct measurements of 
the longitudinal fraction found fo = 0.91±0.39 and 
fo = 0.56 ± 0.31 £ 6j. A recent direct measurement of / + 
set a limit of /+ < 0.18 at the 95% C.L. Q. In addition, 
measurements of the b — > 37 decay rate have indirectly 
limited the V + A contribution in top quark decays to less 
than a few percent ||. However, direct measurements of 
the V + A contribution are still necessary because the 
limit from b — > sj assumes that the electroweak penguin 
contribution is dominant. 

The angular distribution lu of the W boson decay prod- 
ucts with weak isospin ^3 = —1/2 (charged lepton or d, s 
quark) in the rest frame of the W boson can be described 
by introducing the angle 8* with respect to the top quark 
direction 0: 

w(cos6>*) = ^(l-cos 2 r)/ + |(l-cosr) 2 /- 
4 8 



+ -(l + cos6>*) 2 /+ 



(2) 



Due to backgrounds and reconstruction effects, the distri- 
bution of cos 8* we observe differs from uj (cos 8*). How- 
ever, the shape of the measured cos 8* distribution de- 
pends on /_|_ and this dependence can be used to measure 
/+. We do this by selecting a data sample enriched in 
ti events, reconstructing the four vectors of the two top 
quarks and their decay products using a kinematic fit, 



and then calculating cos 8* . This distribution in cos 8* is 
compared with templates for different / + values using a 
binned maximum likelihood method. 

The D0 detector || comprises three main systems: the 
central-tracking system, the calorimeters, and the muon 
system. The central-tracking system is located within a 
2 T solenoidal magnet. The next layer of detection in- 
volves three liquid-argon/uranium calorimeters: a central 
section (CC) covering pseudorapidities \t)\ < 1, and 
two end calorimeters (EC) extending coverage to « 4, 
all housed in separate cryostats. The muon system is lo- 
cated beyond the calorimetry, and consists of a layer of 
tracking detectors and scintillation trigger counters be- 
fore 1.8 T toroids, followed by two more similar layers 
after the toroids. 

This measurement uses a data sample recorded by the 
D0 experiment corresponding to 230 ± 15 pb _1 of pp col- 
lisions at y/s — 1.96 TeV. We consider ti candidate events 
selected in the lepton+jets channel where one of the W 
bosons from t or t decays into an electron or muon and 
a corresponding neutrino and the other W boson decays 
hadronically. The final state is therefore characterized 
by one charged lepton (e or /x), at least four jets (two 
of which are b jets), and significant missing transverse 
energy (# T ). 

Two separate analyses are performed and the results 
are combined. One analysis uses kinematic information 
to select ti events ( "kinematic analysis" ) and the other 
uses b jet identification as well as kinematic information 
in order to improve the signal to background ratio ("6- 
tagged analysis"). A b jet is identified by a displaced 
secondary vertex close to an associated jet [llj . The kine- 
matic analysis vetoes 6-tagged events to simplify the com- 
bination of results with the 6-tagged analysis. In both 
analyses, selected events arise predominantly from three 
sources: ti production, VF+jets production, and multijet 
production where one of the jets is misidentified as a lep- 
ton and spurious $t appears due to mismeasurement of 
the transverse energy in the event. 

The event selection requires an isolated lepton (e 
or n) with transverse momentum px > 20 GeV, no other 
lepton with p T > 15 GeV in the event, Tf, T > 20 GeV, 
and at least four jets. Leptons are categorized in two 
classes, "loose" and "tight," the latter being a subset of 
the first. Loose electrons are required to have \rj\ < 1.1 
and are identified by their energy deposition and isola- 
tion in the calorimeter, their transverse and longitudinal 
shower shapes, and information from the tracking sys- 
tem. For tight identification, a discriminant combining 
the above information must be consistent with the expec- 
tations for a high-pT isolated electron. Loose muons are 
identified using the information from the muon and the 
tracking systems. They are required to have \rj\ < 2.0 
and to be isolated from jets. Tight muons must also 
pass stricter isolation requirements based on the energy 
of calorimeter clusters and tracks around the muon. Only 
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tight leptons are used in the final event selection. Jets 
are required to pass a rapidity 0] cut °f \v\ < 2.5 and, 
in the kinematic analysis, must have px > 20 GeV. The 
requirement that a 6 jet is present significantly reduces 
the background contamination in the 6-tagged analysis 
and allows the use of a lower jet pt cut of pt > 15 GeV 
which increases the efficiency for signal events. 

The top quark and the W boson four-momenta are re- 
constructed using a kinematic fit which is subject to the 
following constraints: two jets must form the invariant 
mass of the W boson, the lepton and the $t together 
with the neutrino p z component must form the invariant 
mass of the W boson, and the masses of the two recon- 
structed top quarks must be equal to 175 GeV. The p z 
component of the neutrino is reconstructed by exploiting 
the fact that the masses of the two top quarks are both 
set to be 175 GeV, and solving the resulting quadratic 
equation for p z 13] . In the case where the two p z solu- 
tions lead to different results of the kinematic fit, the one 
with the lower \ 2 (of the fit) is kept. Among the twelve 
possible jet combinations, the solution with the minimal 
X 2 from the kinematic fit is chosen; Monte Carlo studies 
show this yields the correct solution in about 60% of all 
cases. 

The tt signal events for seven different values of /+, 
/+ = 0.00, . . . , 0.30 in steps of 0.05, are generated with 
the ALPGEN Monte Carlo (MC) program Q for the 
parton-level process (leading order) and PYTHIA for 
simulation of subsequent hadronization. The mass of the 
top quark is set to m t = 175 GeV. As the interference 
term between V — A and V + A is suppressed by the 
small mass of the 6 quark and is therefore negligible , 
these samples can be used to create cos 9* templates for 
any / + value by a linear interpolation of the templates. 
All seven templates from these samples are normalized to 
unit area and a linear fit to the contents of each cos 9* bin 
as a function of / + is performed. This procedure effec- 
tively averages over statistical fluctuations in the gener- 
ated MC samples, thus providing a more precise model of 
the cosO* distribution. The MC samples used to model 
events with W bosons produced in association with jets 
(W+jets) are also generated with ALPGEN, requiring the 
W boson to decay leptonically. The factorization scale 
Q is set to Q 2 = + Y J m T Q. 

To determine the number of multijet background 
events, we compare samples selected with loose and tight 
leptons. Going from loose to tight samples decreases the 
number of events from Ni to N t . The relative selection 
efficiency between the loose and the tight lepton criteria 
is different for true leptons (q) and jets faking an iso- 
lated lepton (ej). We use these efficiencies, known from 
data control samples > to estimate the number of mul- 
tijet background events: N m — (egNg — N t )/(ee — ej). 
The kinematic analysis calculates N m for each bin in the 
cos 9* distribution from the data sample to obtain the 
shape of the multijet cos 9* templates. For the 6-tagged 



analysis, the multijet template is formed from data events 
after the event selection except that the leptons are re- 
quired to satisfy the loose and to fail the tight criteria. 

To discriminate between tt pair pro duction and back- 
ground, a discriminant V is built |l2( using input vari- 
ables which exploit the differences in event topology: Ht 
(defined as the scalar sum of the jet px values), the min- 
imum dijet mass of the jet pairs, the \ 2 from the kine- 
matic fit, the centrality (defined as Ht/He where He 
is the sum of the jet energies) 0, i^mm (defined as 
the distance in rj — 4> space, where 4> is the azimuthal 
angle, between the closest pair of jets multiplied by the 
Pt of the lowest-pT jet in the pair and divided by the 
transverse energy of the reconstructed W boson) [l3j |. 
and aplanarity and sphericity (calculated from the four 
leading jets and the lepton). The last two variables char- 
acterize the event shape and are defined, for example, in 
Ref. 0| . Only the four leading jets in pt are considered 
in computing these variables to reduce the dependence 
on systematic effects from the modeling of soft radiation 
and underlying event processes. All of these variables 
are used for the discriminant in the kinematic analysis. 
Only Ht, centrality, the minimum dijet mass, and % 2 are 
used in the 6-tagged analysis. The discriminant is built 
separately for the kinematic and 6 -tag ged analyses, using 
the method described in Refs. |12lI13i|. The distributions 
of signal (S) and background (B) events in each of the 
above variables are normalized to unity. For each vari- 
able Vi we fit a polynomial to the logarithm of S/B as a 
function of Uj. The discriminant is defined as: 

n , , = expg, MS(vi)/B(vi))ht) , * 

^> V2 ---> exp(E i KS(« i )/BK))]fit) + l ' ( ' 

We select events for which T> > 0.6 in the kinematic 
analysis, and T> > 0.25 in the 6-tagged analysis. These 
values are chosen to minimize the expected statistical 
uncertainty in the measurement of /+ as determined by 
simulations of the analysis. 

We then perform a binned maximum likelihood fit to 
compare the observed T> distribution in the data to the 
sum of the distributions expected from tt, W+jets, and 
multijet events. The number of multijet events is con- 
strained to a Poisson distribution with mean N m . The 
likelihood is then maximized with respect to the number 
of tt, py+jets, and multijet events. We multiply these 
numbers by the efficiency for each type of event to pass 
the V selection to determine the composition of the sam- 
ple used for measuring cos 9* . Table [I] lists the compo- 
sition of each sample as well as the number of observed 
events in the data. The cos 9* distribution obtained in 
data after the full selection is shown in Fig. ^ f° r the 
kinematic and in Fig. for the 6-tagged analysis. 

A binned maximum likelihood fit of signal and back- 
ground cos 9* templates to the data was used to measure 
/+. We compute the binned Poisson likelihood (L(/+)) 
of the data to be consistent with the sum of signal and 
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TABLE I: Number of events observed for each component 
(signal and backgrounds) of the kinematic and 6-tagged sam- 
ples, and the number of data events, after the cut on the 
discriminant T> discussed in the text. The fits are made be- 
fore the final cut on T>, so the sum of the components need 
not agree exactly with the observed numbers of events. 



TABLE II: Systematic uncertainties on /+ for the two inde- 
pendent analyses and for the combination. 



Event Class 


Kinematic 


b- tagged 


tt 


16.5 ±5.8 


40.8 ±8.1 


W+jets 


14.3 ±3.0 


11.5 ±4.1 


Multijet 


5.0 ±2.1 


1.5 ±0.5 


Data 


35 


52 



C 15 

> 
LU 

10k 



-•- D0 data 
I I tt-> l+jets (V-A) 
- - - tt-> l+jets (V+A) 
background 




FIG. 1: cos 8* distribution observed in the kinematic analy- 
sis. The standard model prediction is shown as the solid line, 
while a model with a pure V + A interaction would result in 
the distribution given by the dashed line. 



background templates, normalized to the numbers given 
in Table at each of the seven chosen / + values. In 
both analyses, a parabola is fit to the — ln[L(/+)] points 
to determine the likelihood as a function of / + . 

Systematic uncertainties are evaluated in ensemble 
tests by varying the parameters (see Table [HJ which can 
affect the shape of the cosf?* distributions or the rela- 
tive contribution from the three sources (ti, VF±jets and 
QCD). Ensembles are formed by drawing events from a 
model with the parameter under study varied. These are 
compared to the standard cos 0* templates in a maximum 



£ 20 
c 

■2 15 



-•- D0 data 

I I tt^ l+jets (V-A) 

tt-> l+jets (V+A) 

[~x| background 




0.5 1 

cose* 



FIG. 2: cos 8* distribution observed in the 6-tagged analy- 
sis. The standard model prediction is shown as the solid line, 
while a model with a pure V + A interaction would result in 
the distribution given by the dashed line. 



Source 


Kinematic 


6-tagged 


Combined 


Jet energy calibration 


0.03 


0.04 


0.04 


Top quark mass 


0.04 


0.04 


0.04 


Template statistics 


0.05 


0.02 


0.03 


6-tag 


0.03 


0.02 


0.02 


ti model 


0.01 


0.02 


0.02 


W±jets model 


0.01 


0.01 


0.01 


Sample composition 




0.02 


0.01 


Calibration 


0.01 


0.01 


0.01 


Total 


0.08 


0.07 


0.07 



likelihood fit. The average shift in the resulting / + value 
is taken as the systematic uncertainty and is shown in 
Table ITT1 The total systematic uncertainty is then taken 
into account in the likelihood by convoluting the latter 
with a Gaussian with a width that corresponds to the 
total systematic uncertainty. 

The top quark mass and the jet energy calibration 
(JEC) are the leading sources of systematic uncertainty. 
The mass of the top quark has been varied by ±5 GeV 
with respect to m t = 175 GeV and the JEC by ±ler 
around the nominal value. The statistical uncertainty 
on the cos 6** templates has been taken as a systematic 
uncertainty. It is estimated by fluctuating them accord- 
ing to their statistical uncertainty. Uncertainties in the 
modeling of the 6-tag algorithm lead to uncertainties in 
the flavor composition of the W+jets background and 
in the cos 0* distribution itself due to the pr and rj de- 
pendence of the 6-tag algorithm . An uncertainty in 
the flavor composition translates into a different shape 
of the cos 0* distribution and a difference in the signal 
to background ratio. In order to estimate the systematic 
uncertainty due to gluon radiation in ti events, an alter- 
native signal sample of t£+jet has been generated with 
ALPGEN, and mixed with the default ti sample using the 
leading order cross sections for both processes. Effects 
of the choice of factorization scale Q in the generation 
of the IU±jets events have been evaluated by using a 
sample where Q 2 = (pt) 2 Q- There is a systematic un- 
certainty due to the final sample composition obtained 
by the fit to the discriminant T>. The kinematic analysis 
treats this uncertainty as a statistical uncertainty and in- 
cludes it in the definition of the likelihood as described in 
Ref. 0| while in the 6-tagged analysis this uncertainty is 
studied by changing the compositions within their errors. 
The difference found between the input /+ value and the 
reconstructed /+ value in ensemble tests is taken as sys- 
tematic uncertainty on the calibration of the analysis. 

The result of the maximum likelihood fit to the cos 9* 
distribution observed in the data is shown in Figs. OJa) 
and (b) for the kinematic and 6-tagged samples, respec- 
tively. The statistical uncertainties from the two indi- 
vidual analyses are 0.22 for the kinematic and 0.17 for 
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FIG. 3: — InL curve obtained in the a) kinematic analysis, b) 6-tagged analysis, and c) kinematic and 6-tagged analyses 
combined. The dashed line includes only the statistical uncertainty while the solid line also includes the systematic uncertainties. 
The physically allowed region for /+ is indicated by the grey area. 



the 6-tagged analysis. The — ln[L(/ + )] curves for the 
kinematic and ^-tagged measurements are combined, as 
shown in Fig. EJc). The systematic uncertainties are as- 
sumed to be fully correlated except for the systematics on 
calibration of the individual analyses which are uncorre- 
cted, and the Monte Carlo model systematics which are 
partially correlated. Assuming a fixed value of 0.7 for /o, 
the combined result for / + is: 

/+ = 0.00 ± 0.13(stat) ± 0.07(syst). (4) 

The observed combined statistical uncertainty (0.13) is in 
good agreement with the expectation (0.12) inferred from 
ensemble tests. We also calculate a Bayesian confidence 
interval (using a flat prior distribution which is non-zero 
only in the physically allowed region of / + = 0.0 — 0.3) 
which yields 

/+ < 0.25 @ 95% C.L. (5) 

The W boson positive helicity fraction /+ that we 
have measured in ti decays in the lepton+jets chan- 
nel is consistent with the standard model prediction of 
/+ = 3.6 x 10~ 4 0. 
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